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Abstract 

We develop a series expansion of the plasma screening length away from the classical limit in 
powers of h'^. It is shown that the leading order quantum correction increases the screening length 
in solar conditions by approximately 2% while it decreases the fusion rate by approximately 0.34%. 
We also calculate the next higher order quantum correction which turns out to be approximately 
0.05%. 
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I. INTRODUCTION 



Salpeter^ wrote a seminal paper more than half a century ago, concerning screening effects 
on thermonuclear reaction rates. He made the basic point that screening effects are small 
at the center of the sun. 

There has been renewed interest more recently in utilizing the sun as a source of neutrinos 
to test the standard model of unification of electro-weak forces. The measured neutrino flux 
deviates from predictions of the standard model by a factor of two^. The measurement un- 
certainty is^'^ ~ 1%. Therefore it would be meaningful to quantify the theoretical estimate 
with equal precision. The structure and dynamics of the sun are complex^. Various phenom- 
ena need to be identified and estimated correctly. Screening of Coulomb repulsion between 
nuclei at extremely short distances is one of them. Many calculations of screening have 
been made since Salpeter's original paper, attempting to refine the degree of screening^'^'^, 
dynamic effects^, quantum fluctuations^'^, etc^°. Here we shall focus on quantum correc- 
tions to screening. The most sophisticated calculation of this effect is that of Gruzinov 
and BahcalP. In this paper, the electronic density matrix was evaluated accurately using 
Feynman's formulation in terms of a Schroedinger equation with the inverse temperature 
playing the role of imaginary time. Fermion statistics are ignored due to the high solar tem- 
perature. They sustain Salpeter's original conclusion that quantum corrections are minor. 
These calculations are essentially correct, but cannot estimate in a systematic fashion the 
next higher order quantum correction. We shall correct that deficiency in this paper. 

This paper was written for the sake of completeness, since the super-Kamiokande exper- 
iment has been successful in obtaining evidence for neutrino mass (see Hosaka et al^^ for 
recent results). Nevertheless, the results of this paper may yet prove useful for more precise 
quantitative interpretation of stellar experimental data^^. 

II. NON-RELATIVISTIC QUANTUM CORRECTIONS 

We shall treat ions in a thermonuclear plasma as classical objects, while applying a quan- 
tum treatment to electrons. The resulting partition function is evaluated as a deviation from 
the classical limit for electrons. A series expansion is developed in powers of a dimensionless 
ratio involving fi^. 



Let us begin with the classical Poisson-Boltzmann equation for a single species of ions 
and electrons: 



- = 47rp 

P = P+ + P- 
p+ — e n Z ex^{—Ze(f)/kBT) 

p_ ^-enZ ex-p{e(f)/kBT) (1) 

where e is the magnitude of the electronic charge, is Boltzmann's constant, n is the 
average number density, Ze is the ionic charge, and T is the temperature of the system. We 
shall work in the linear regime, which is expected to apply in the solar interior^, by retaining 
only terms first order in 0: 
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where Aq is the classical screening length. 

The quantum-mechanical version of this approximate Poisson-Boltzmann equation for a 
single species of ions and electrons may be written in analogy with Eqn. 1: 



- = 47rp 

P ^ P+ + P- 

p_ =-e|^({f})r (3) 

where ifj is the many-body quantum wave-function for electrons, and {f} refers collectively 
to the electrons in the system, and (p is the electrostatic potential. 

We now invoke the following scaled variables, in order to ease subsequent calculations: 



' =A 

Note that the first of Eqns.4 shows that we are using UbT as the energy scale. The electro- 
static potential is then given by: 

V''0 = (0 + 47rr|^|2-Z-^) (5) 
This equation may be obtained from a Lagrangian density: 

>Co ^-\\V4>\^-v{4>.^) 
v{4>,^) =^4>' + ^7i4>m'-Z''^ (6) 

The corresponding Hamiltonian density can be easily derived: 

^o = ^|V0|^ + ^(0,^) (7) 
We will now introduce second-quantized notation to deal with the statistics of electrons: 

v{^, i>) ^ v{^, vi±) = - z-'4> + 47r0r(viiv^+ + ipli'-) (8) 

where iIj± are Grassmann variables, and the subscripts refer to the spin of the electrons. 
The co-existence of Grassmann variables and scalars in Eqn.8 is not problematic, since we 
shall use this discussion solely to define a partition function for the entire system. And soon 
thereafter we shall integrate over the electron degrees of freedom, so that only a functional 
involving the scalar potential survives. 

The total Hamiltonian Ti for the system, including the quantum-mechanical part for the 
electrons is: 



hq ^Agdw+i'+iw-r) 



(9) 



The quantum correction has been encapsulated in the following dimensionless parameter: 



where m is the mass of the electron. 

Since solar temperatures are ~ 0{lkeV), and the rest energy of the electron is 0.55Mey, 
it follows that the non-relativistic approximation employed in Eqn.9 is valid. 

The partition function may be written in scaled variables as: 



where it is understood that ksT = 1 in the units we are using. 

Note that the total number of electrons is associated with each ion of charge Ze is Z, 



expectation value). We will indicate shortly how one may impose this constraint upon 
the system using a Lagrange multiplier. In condensed matter physics, this is done via the 
electronic chemical potential. In our problem, it will turn out to be more convenient to 
institute this constraint via a functional involving just the electrostatic potential. 

Note that the parameter F is analogous to the usual plasma parameter. It is much 
less than one for solar conditions. For solar conditions, viz., a density of 100 g cm~^, 
T = 15 X lO^K, Z = 2, it turns out that Aq ^ 0.28ll, F ^ 0.04. The value of Aq turns out 
to be approximately 0.032. This is already an indication that quantum corrections will be 
small. In making these estimates, we have assumed a Helium plasma with thermodynamic 
properties similar to those at the center of the sun^. 

The quadratic nature of the energy functional in Eqn.ll allows us to perform the func- 
tional integration over the Grassmann variables associated with the electronic degrees of 
freedom^^, allowing us to obtain: 




(10) 




(11) 



and is obtained via < {ipl^tp^ + ip'Lip^) >= Zn (where the angular brackets indicate an 



Z j exp(- j ciV((l/2)|V0|2 + (1/2)02 _ ^-i^)) Det(JP) 
Det(:r) =exp(Trln(:r)) 

JT = -Aq + 47rr0 (12) 

Having integrated over the electronic degrees of freedom, we are left with an effective 
energy density in terms of the electrostatic potential alone. We choose to impose charge 
neutrality, which was discussed just below Eqn.ll, via a Lagrange multiplier by making the 
following addition (A7i) to the energy density: 

AH = 47ri/0 (13) 

Here u may be interpreted physically as a uniform charge density, which will be adjusted 
to ensure the overall charge neutrality. 

We need to evaluate the determinant of the operator obtained in the process of of perform- 
ing the quadratic functional integral over fermionic variables. This is conveniently performed 
in Fourier space: 

^ HATrVm + Aq e) (14) 

where 0(A;) is the Fourier transform of 0. 

Now the estimates below Eqn.ll indicate that 47rr >> Aq near the center of the sun, so 
we propose a series expansion in powers of Aq: 

HAT^Vm + Aq k') ^ HAT^Vm) + -^^i;-, - i^iir^ + - (15) 

47rr0(/c) \ATiV(t){k) j 

Furthermore, we will conform to the linear screening limit which is expected to apply in 
the solar interior^. By this we mean we shall seek an expansion of the determinant obtained 
above, to the quadratic order of the scalar potential.- higher order terms can be accounted 
for using diagrammatic techniques. We have the freedom to choose the value of the potential 
around which to perform the expansion. We shall choose this value to be 0o, such that in 
the limit that ^ ^ 0, we recover the standard expression for the screening length obtained 
from the linearized Poisson-Boltzmann equation. The power series to quadratic order yields: 



ln(47rr0(A;) + Aq k^) ^ ln(47rr) 

+(ln(0o) + im - </'o)/0o - - 0o)V(20^)) 

+^ (1 - im - 0o)/0o + {m - 0o)V0o) (16) 

The constant terms are unimportant as they can be absorbed into the normahzation 
constant. To be consistent, the rest of the effective energy density must also be expanded 
around 0o- The net resuh would yield, in addition to a quadratic term, a term linear in <l>. 
The coefficient of this linear term is an effective background charge density, which must be 
zero in our neutral system. The coefficient can be set to zero by adjusting appropriately 
the Lagrange multiplier, defined in Eqn.l3. 

With the proper charge neutrality constraint imposed, the screening length expression 
can be matched with the linearized classical Poisson-Boltzmann Eqn.2 when — * 0, by 
setting 00 = V^- Then the energy density may be written to leading order in Aq in Fourier 
space as follows: 

1 9 A 1 

'H^-k'{l + ^Mk? + + l/Z)m' (17) 

The corresponding Lagrangian density in real space is: 

1 OA 1 

>C = -2 (1 + - + V^)^(r1^ (18) 

The equation of motion then follows: 



y2^ _^^-2^ 




2Aq 



In dimensionful units, the corrected screening length is: 



- V 47r(Z^ + Z)ne- 

Note that without quantum corrections, the screening length is about 0.28lA. Quantum 
fluctuations increase the screening length by ~ 2% at solar conditions, deflned earher. The 



classical enhancement factor of the fusion rate is 1.17. It is reduced shghtly to 1.16 via 
quantum corrections. The next higher order quantum correction to the rate from our theory 
(expanding Eqn.l6 to 0{Aq)) turns out to be approximately 0.05%. 

The numerical values obtained for the screening length and the Salpeter rate enhancement 
factor have been encapsulated for solar conditions in Table I: 

Quantity Quantum Correction 

Screening Length 2% 
Rate enhancement 0.34% 

TABLE I: Comparison of classical and quantum-corrected quantities in solar conditions. 

Hence the leading order quantum correction decreases the fusion rate by about 0.34% for 
the conditions chosen. The numerical estimates provided above are in fair agreement with 
those found in the literature^ 

III. CONCLUSION 

Systematic quantum corrections to screening in thermonuclear fusion were derived in 
powers of and estimated for solar conditions. Leading order corrections were shown to 
be less than 0.34% under solar conditions, while the next leading order term is ~ 0.05%. 

Our corrections are consistent with those previously obtained by Gruzinov and BahcalF. 
They complement the results of Brown et al^ who show that classical, non-linear effects on 
screening are small. 
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